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Starting from the Wigner-Moyal equation coupled to Poisson’s equation, a simplified set of equa¬ 
tions describing nonlinear Landau damping of Langmuir waves is derived. This system is studied 
numerically, with a particular focus on the transition from the classical to the quantum regime. In 
the quantum regime several new features are found. This includes a quantum modified bounce fre¬ 
quency, and the discovery that bounce-like amplitude oscillations can take place even in the absence 
of trapped particles. The implications of our results are discussed. 
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I. INTRODUCTION 

Numerous aspects of quantum plasmas have been investigated during the last decade, see e.g. mm- Basic features 
such as electron degeneracy and particle dispersive properties have been studied in some detail m Exchange effects 
mm, the magnetic dipole force and other contributions from the electron spin mm have also been examined, 
including relativistic effects mm- Systems of interest in this context include e.g. quantum wells m , laser-plama 
interaction on solid density targets m, and astrophysical plasmas M- Works on quantum plasmas also have relevance 
for recent applications in plasmonics 19 and spintronics [20] . 

In the present paper we will consider the influence of quantum effects on the nonlinear regime of Landau damping, 
based on the Wigner-Moyal equation [2] that accounts for particle dispersive effects. Previous works in this area 
mm have deduced that quantum effects can suppress the nonlinear bounce oscillations of Langmuir waves and 
turn the evolution into basic linear damping. Moreover, quantum corrections to electron holes in phase space that 
may form as a result of wave-particle interaction have been calculated in Ref. m , and the quasilinear theory of the 
Wigner-Poisson system has been studied in Ref. [24]. However, many of the details of nonlinear Landau damping 
have not been studied before in the quantum regime. Making analytical approximations applicable for a resonance in 
the tail of the distribution, we first simplify the Wigner-Moyal equation coupled to Poisson’s equation into a system 
that is more easy to solve numerically. This system is shown to fulfill an energy conservation law, and reduces to a 
previously studied system [25] in the classical limit. A systematic study of the transition from classical to quantum 
behavior then reveals several new features. This include a quantum modification of the bounce frequency, a new 
condition for the quantum suppression of the nonlinear regime and the discovery that bounce-like oscillations can 
take place even in the absence of trapped particles. 

Similar to most cases, the conditions needed for quantum effects to be important in our problem include a high 
plasma density and a modest plasma temperature. However, while the scaling with temperature and density is the 
same as usual mm, the precise numerical values are to some extent relaxed compared to the standard expressions in 
case the resonance lies in the tail of the distribution. The general conclusion is that the properties of wave-particle 
interaction are more easily influenced by quantum effects than the ordinary fluid properties. A concrete example is 
provided in the final section, and the significance of our results are discussed. 


II. BASIC EQUATIONS AND DERIVATIONS 

Our starting point is the Wigner-Moyal equation [2], which reads 
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Here / is the Wigner-function, is the electrostatic potential, q (=—e) is the electron charge, m is the electron mass 
and h = 2irh is Planck’s constant. Eq. 0 is combined with Poisson’s equation 

- v 2 $ = — [ fd 3 v (2) 


to give us a closed set. Eq. 0 applies for electrostatic fields and do not account for the spin of the electrons or 
exchange effects. For generalizations to electromagnetic fields and spin effects, see e.g. Ref. [123 an d f° r generalizations 
to include exchange effects, see e.g. m- Here we focus on the problem of Langmuir waves in an unmagnetized plasma, 
in which case the omissions of electromagnetic effects and spin effects are trivially justified. Furthermore, we will 
consider the case of a moderate plasma density with huo p It is then safe to exclude exchange effects m ■ In 

fact, the condition hw p ksT is often used also to disregard the particle dispersive quantum effects included in Eq. 
0 . Such a condition is indeed an appropriate one for neglecting e.g. the particle dispersive contribution to the real 
part of the Langmuir dispersion relation and for neglecting the quantum contribution of Eq. 0 in many other cases 
mm ■ However, as we will see below, the quantum contribution to wave particle interaction can be crucial even in the 
regime fouj p <C A useful starting point to see this is to rewrite Q in the alternative way [26] 
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Here the arrows on the operators indicate in which direction they are acting, and the sinus-operator is defined in 
terms of its Taylor-expansion. Let us first consider linear plane wave solutions oc exp i(kz — c ot), and estimate the 
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relative importance of quantum corrections. Letting / = / 0 + fi exp i(kz — ujt) we can use the classical linear solution 


A- ( 4 ) 

m{uj - k z v z ) dv z 

as a means to illustrative the importance of quantum corrections. Comparing the magnitude of the first order quantum 
correction in the Taylor expansion with that of the classical term, we note that in the bulk of the velocity distribution 
(i.e. for v ~ v t ) the quantum corrections are small provided Qbuik 1, where 
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However, looking at velocities close to the wave-particle resonance we must let uj — k z v z be small, of the order 
uj — k z v z rsj where yz, is the linear damping rate. Close to the resonance the relative importance of the quantum 
terms are given by the parameter Q res , which is 
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For a resonance in the tail of the distribution y^/fc v t m which case we may have Q res ^ 1 at the same time as 
Qbuik 1. This regime will be the focus in what follows. As a result it suffices to solve the Vlasov limit of Eq. © 
in most of velocity space, but close to the resonance we need to solve the full Wigner equation. A similar approach 


applies for the small amplitude approximations. Provided qk&/mnujv t < 1 we have df\/dv z <C \dfo/dv z \ unless we 


are close to the resonance. We will assume qk&/mujv t <C 1 to hold, and thus for most of velocity space we can solve 
the linearized Vlasov relation. However, close to the resonance where we can have nonlinear wave-particle interaction 
we must solve the Wigner equation without making linear approximations. 

For a resonance in the tail of the distribution the number of particles contributing to the nonlinear interaction are 
relatively few. As a result there will be only minor harmonic generation of the electric field. Consequently we will 
use the ansatz of a slowly time-varying plane wave <f> = <!>(£) exp i(kz — cut) + c.c. for the potential, where c.c. stands 
for complex conjugate In general the Wigner function will be a periodic function, that will be represented as 


/ = /o(v) +6fo(v,t) + 
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Next we divide the velocity space into the resonant region [v z — Sv re s, v z + £u res ] and the non-resonant region containing 
the complementary part. When evaluating the charge density in O) we thus use f fd 3 v = f nr fd 3 v + f res fd 3 v where 
the subscripts nr and res denotes the nonresonant and resonant regions, respectively. Since the linear Vlasov equation 
applies in the nonresonant region we can use 
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Noting that the second term on the right hand side of 0 is a small correction, proportional to the slowly varying 
amplitude, we can use the lowest order approximation (i.e. dropping time derivatives on the amplitude) to convert it 
to a term proportional to d&{t)/dt. Combining |?| with ([i]) we then deduce 


d$(t) 

dt 


Q 

£ok 2 dD/duj 



( 8 ) 


where D(u;, k) = l+(q 2 /kmso) f nr (dfo/dv z )d 3 v/(uj—k z v z ). Note that the dispersion function D(uj, k) has a component 
of arbitrariness in the definition, since it depend on the width of the resonant region £u res . Nevertheless in the 
derivation of (|8| we have taken D(u,k) = 0 to hold by definition (i.e. w(k) is defined by this relation), and in case 
a possible frequency shift occurs due to this, it is included in the time-dependence of . It should be stressed that 
although a large number, of harmonics / n (v,£) might be needed to solve for the Wigner function in the resonant 
region, it is only the first harmonic /i that contributes in (|8|. 
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Next we need to solve for the Wigner function in the resonant region. We restrict ourselves to a Maxwellian 
background distribution Fq(v), and introduce the lD-Wigner function 


g — Gq(v z ) + go(v z ,t) + 
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where Fq(v) = Gq{v z ) exp[(—and / n (v,£) = g n {y z , t) exp[(—— r^)/?^]. This ansatz is then substituted 
into ([3|. Since the gradient operator becomes ±ikz (as the spatial dependence of the potential is exp (±ikz)) a useful 
formula is 

2m . ( ihk d\ ±im[g n (v z + Hk/2m) - g n (v z - hh/2m)\ , 

) = -B-' (10) 

With the help of this formula, and subsituting the ansatz © into ©> the following set of coupled equations are 
deduced 
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where Suj(v z ) = u; — kv z and the star denotes complex conjugate. Using / reg /id 3 ?; = ttv 2 f res gidv z in Q we see that 
Eqs. © and © constitute a closed set. The equations agree with those studied by Ref. [25] in the classical limit 
h —>> 0, provided the collisional frequency is put to zero in that work. It is clear that the quantum features are encoded 
in the velocity shift hk/2m. An important thing to note is that we can use Gq{v z + hk/2m) — Go(v z — hk/2m) ~ 
[dGo(v z )/dv z ]hk/m for hw/mv 2 1, whereas similar approximations are not applicable for the perturbed Wigner 
function, as the quantities g n varies on a much shorter scale length. 

Before we can proceed we must put constraints on the parameter Sv res . For the resonant region to cover a sufficient 
amount of resonant particles in the linear regime we need the condition d?; res 7 B /k. Moreover, to cover particles 
that are closed to be trapped in the potential well with a sufficiently large margin we need d?; res uos/k, where 
uj b = (qk 2 <&/m) 1 / 2 is the bounce frequency of trapped particles (with this choice the the resonant region is at least 
an order of magnitude larger than the region of trapped particles). Finally, to cover the quantum effects properly we 
need Sv res hk/2m. At the same time the calculation scheme is based on the resonance region being much smaller 
than the thermal velocity, and hence we need d?; r es <C v t . If we sharpen this condition slightly and limit ourselves to 
Sv res kv 2 /u we may take [dGo(v z )/dv z ] as constant in the resonance region (= [dGo(v z )/dv z \ UJ /k)i which simplify 
some of the technical aspects. For a resonance that lies in the tail of the distribution we can have kv t /^L ~ 100, 
in which case it is easy to fulfill all conditions simultaneously, unless the nonlinearity or the quantum effects are 
extremely strong. Importantly, the numerical solutions presented below are not dependent on the precise choice of 
Sv res , as long as the above conditions are fulfilled. 

Next we introduce normalized variables. Choosing normalized time as 7 ^, normalized velocity as kv z /^L, normal¬ 
ized potential as qk 2 $/^\m and normalized harmonics of the Wigner function as 

kg n 
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the coupled equations become 
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and 
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where n > 2 in Eq. (16) and the quantum velocity shift is 5v q = hk 2 /2m'yL- We have omitted indices denoting 
normalized variables for notational convenience. Provided the conditions for the resonance region presented above 
is fulfilled, we may neglect perturbations at the boundary, i.e. use the approximation g n (±Sv res ) « 0 where n = 
0,1,2....This property holds in the vicinity of the boundary (where the vicinity means a velocity of the order 5v q ), 
which means that 


/ 9 o(v + 6v q )gZ(v)dv « / go(v)g{{v - dv q )dz 

J res J res 
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When Sv res is chosen large enough such that Eq. (18) is fulfilled, as well as similar types of approximations involving 
# n , the system posses an energy conservation law 
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where VE to t is the total energy. The first term of Eq. (19) represents the wave energy (including the kinetic energy in 


the non-resonant region), whereas J reg | g n \ dv represents the particle energy in the resonant region of each harmonic. 
Note that the zeroith harmonic get an extra factor 1/2, which is related to go being real. In the numerical calculation 
made in the next section Eq. (fl9|) has been used as test of the numerical scheme and a confirmation that Sv res has been 


chosen large enough. In particular it should be stressed that the quantities J reg | g n \ dv is not sensitive to the exact 
choice of Sv res , since the integrand falls of rapidly away from the resonance. This can tested by varying the parameter 
Sv res in the numerical code. We find that the relative values of the wave energy and resonant energy changes around 
0.001 — 0.002 when the value of Sv res is changed a factor 1 — 3 within the bounds of the strong inequalities. 

An advantage with the above system ( pT| )- is that solving the equations numerically we can follow the evolution 
taking time-steps that are larger than the inverse plasma frequency, as the equations contain only the slow time-scales, 
as opposed to the original Wigner-Moyal equation ©• Moreover, we only need to solve the equations in a small part 
of the velocity space, close to the resonance, which also makes Eqs. easier to solve numerically. Finally, the 

spatial dependence is solved for analytically in (fl4l)-([rF|), which also simplifies the numerics. 


III. REMARKS ON THE LINEAR THEORY 


Before we start with a numerical study, let us first make a few comments regarding the linear theory. As the 
equation stands, making a linearization of Eqs. ([l4|) - (p~7|) completely removes the quantity Sv q which encodes the 
quantum properties. Eq. (15) can then be integrated according to 
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For certain initial conditions gi(t = 0) the integrals in Eqs. (14) and (20) can be solved analytically (see e.g. Ref. 
m for details), and for these cases we indeed have linear damping 4> oc e 1 which is written e lLt in non-normalized 
units, where 
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The reason no effects due to the quantum treatment is seen here is the assumption of a modest quantum regime 
kv 2 /uj Sv res hk/2m. These conditions apply for a resonance in the tail of the distribution, when the wavelength 
is not too short. Whenever these inequalities hold we can also use 


Gq(v z + hk/2m ) — Gq{v z — hk/2m) 
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which is the reason the finite difference is replaced by the classical expression containg a derivative in Eq. (21). 
Since go, gi,g 2 ••• varies on a much shorter scale length in velocity space as compared to Go, we note that a similar 
approximation cannot be applied for these quantities. We stress that in a general scenario (e.g. for a beam-plasma 
system or for the regime hk 2 /m ~ uj) the approximation (J22| may give an incorrect description of linear Landau 
damping, but in our case the linear regime is classical to a good approximation. It shold also be noted that the 
evolution of g\ exhibits phase mixing due to the factors e~ lvt (see also Fig. 1 of the next section) whether or not the 
linear damping expression contains a finite difference or a derivative. Hence the quantum influence on the regime of 
linear Landau damping tend to be relatively modest also when the approximation (22) is avoided. 


IV. NUMERICAL SOLUTIONS 




v 


Figure 1. The real part of part of g\ evaluated in the small amplitude regime for t = 2 (upper panel) and t — 8 (lower panel). 
The term ivg\ of Eq. (15) leads to phase mixing, i.e. the development of increasingly small scales in velocity space. It should 
be noted that quantum effects does not counteract this behavior for small amplitudes. 



Figure 2. A qualitative illustration of the trapped particles. For a fixed wave amplitude particles trapped particles in the wave 
field have discrete energies, with an energy step of the order of huts- When this becomes comparable to only a few states 
are trapped, and quantum mechanical effects become important to describe the evolution. 


The system @-([17]) has been studied numerically using a staggered leapfrog finite differencing technique. Two 
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classes of initial conditions have been used. In the first case gi(v, t = 0) = 4>(t = 0)/(iv +1) and go(t = 0) = g n (t = 0) 
= 0 (for n > 2). This has the advantage that the Wigner function approaches the value outside the resonance reegion 
for v 1. The other choice is to put also gi(v,t = 0) = 0 such that the initial Wigner function is identically zero 
in the resonant region. It turns out that the evolution is almost completely independent of this difference in initial 
conditions, if the initial amplitude 4>(t = 0) is increased in the latter case such as to make the initial energy equal in 
the two cases. In what follows all numerical results will refer to the case with gi(v, t = 0) = 4>(t = 0 )/{iv + 1). Picking 
a small initial amplitude = 0) < 1 the system shows a damping 4> oc e~ l as expected due to the normalization 
of the time variable. As noted in the previous section the evolution is independent of the quantum parameter 5v q , 
which is related to the assumption fruj/mv 2 <C 1 that was made in the derivation. Fig. 1 compares the real part 
of g\ for different times. The evolution shows the evolution towards smaller scale lengths in velocity space due to 
phase mixing. When nonlinearities come into play this process plays part in increasing the relative importance of 
quantum effects in wave-particle interaction, as mathematically the system ( p~4| )- © differs from the classical limit 
once the scale length of g n in velocity space is smaller than Sv q . However, a more physical way to understand why 
quantum effects become important more easily in the nonlinear regime than in the linear regime is shown in Fig. 2. 
Here as a rough approximation we have described the discrete eigenstates of trapped particles as that of an harmonic 
oscillator where the eigenfrequency is the bounce frequency ujb- Naturally this is rather crude, as the electrostatic 
potential is only harmonic for the lowest energy states (if at all), and moreover the potential can vary dynamically to 
a smaller or larger degree depending on the parameters of the problem. Still the simple picture in Fig 2 is sufficient 
to identify one of the key parameters R tr , which is the ratio of the trapping potential over the energy quanta of 
trapped particles, R tr = q&/hujB- If the initial values have R tr 1 the evolution of the wave amplitude 4>(£) is 
classical to a good approximation. On the other hand, when R tr decreases towards unity the discrete energy states 
of the trapped particles will modify the evolution significantly. In terms of our normalized variables we note that the 
quantum trapping parameter can be written as R tT = q^/hujB = muiB/hk 2 = 4* 1 / 2 /Sv q . For R tr < 1 we will not have 
trapped particles. However, as we will see below the close coupling between the existence of trapped particles and 
nonlinear evolution that holds in the classical regime does not generally apply in the quantum regime. 



Figure 3. The evolution of 4>(t) for initial amplitude 4>(t = 0) = 4 and 5v q = 0.5, 4, 8 and 20. The value Sv q = 0.5 gives a 
classical evolution to a good approximation, whereas for 5v q = 20 the nonlinearities are almost completely suppressed by the 
quantum effects. 


Next we investigate how the value of the quantum parameter 5v q affects the evolution of 4>(£) in the nonlinear 
regime. Keeping the normalized initial amplitude equal to 4>(t = 0) =4 (which means that the (initial) bounce 
frequency is ujb = 2y£,), we follow the evolution up to t = 17.5 for various values of Sv q . The result is shown in Fig. 
3. For Sv q < 1 the evolution of the wave amplitude more or less coincides with the classical case. This means that 
an initial drop in amplitude is followed by oscillations with a frequency of the order of the bounce frequency ujb • As 
ujb is not a constant, it is no surprise that these oscillations are slightly irregular. This result is in agreement with 
previous works on the classical case, see e.g. [2> , I23 EH] . When the quantum parameter 6v q is increased, more of 
the wave energy is converted to the particles (the initial amplitude drop is larger) and the period of the amplitude 
oscillations become longer. Eventually when Sv q = 20 the initial amplitude drop is so large such that the evolution 
has become almost completely linear. Keeping Sv q = 20 and instead varying the initial amplitude, the evolution 4>(t) 
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Figure 4. The evolution of 4>(t) for Sv q = 20. The inital amplitudes are <t(£ = 0) = 4, 12, 20 and 40. 

for various values of <F(£ = 0) is shown in Fig 4. The same qualitative features as in Fig 3 can be seen. That is 
the amplitude oscillation gets a lower frequency with decreasing $(t = 0), and the initial drop in wave amplitude 
is larger, until eventually for small enough initial amplitude the nonlinearities are suppressed and we obtain linear 
damping. A more quantitative analysis based on Figs 3 and 4 reveals that amplitude oscillations have a frequency 
Wamp ~ (u%/(l + Sv q ) — 7 I/ 2 ) 1 / 2 Classically Sv q = hk 2 /2m'yL is small and uj amp ^ cj b for ujb/il 1. When 
the quantum velocity shift Sv q becomes comparable to the characteristic scale of g\ the frequency of the amplitude 
oscillations decreases. Still the amplitude may oscillate nonlinearly even for large Sv q and R tr <C 1 (in the regime of 
no trapped particles). The transition to linear evolution occur when the expression for u; am p becomes begative, in 
which case the amplitude oscillation frequency decreases too fast for the nonlinear oscillations to get started. 

Next we look closer on some of the details of the numerical results. As shown in Figs 3 and 4, the evolution of the 
wave amplitude changes rather smoothly with 5v q . However, the various harmonics of the Wigner function g n is more 
sensitive to the change in Sv q . The discrete structure involving momentum changes hk = 2 mSv q can be seen directly 
in the Wigner function, perhaps most clearly in g$. In a very rough sense there is a relation between the energy 
states displayed in Fig. 2 and the momentum structure shown in the upper panel of Fig. 5. As shown in Fig 5 where 
go(v, t = 10) is plotted in the classical and quantum mechanical regime, respectively, there is a very distinct difference 
between the quantum mechanical and classical regime. For the former case clear dips for v ~ —2 Sv q and v ~ — 5v q 
can be seen, as well as peaks for v ~ Sv q and v ~ 2 Sv q . It should be noted, however, that a discrete structure in the 
momentum dependence can be seen even in the absence of trapped particles. 

Another important aspect is the convergence in the sum over harmonics g n , which is much faster in the quantum 
regime. A comparison of the classical and quantum regime is displayed in Fig. 6. We see that the relative amount 
of energy in the third harmonic f |^ 3 (t)| 2 dv/W tot is changed dramatically when Sv q is changed from 0.5 to 4. The 
convergence in the sum over the harmonics is determined by the trapping parameter R tr . In the classical and nonlinear 
regime with R tr > 1 we need to include harmonics up to #io, or even more. Besides the value of R tr the number of 
harmonics needed depends on the degree of the nonlinearity and how long the evolution is followed. In the quantum 
regime when R tr > 1 including up to g% is typically more than enough to get convergence. This applies even when 
Wamp 7 l and the evolution is strongly nonlinear. While the trapping parameter R tr determines much of the 

properties of the Wigner function in the resonant region, it is the parameter R n \ = + Sv q )] = <f>/(l + Sv q ) 

that determines the degree of nonlinearity. In the nonlinear regime R n \ >1/2 this parameter determines the frequency 
of the nonlinear amplitude oscillations, and we have cj 2 mp = (-Rni — l/2)y£- For small R n \ (in practice R n \ < 0.25) the 
evolution is linear to a very good approximation. The reason why R n \ determines the nonlinearity can be understood 
roughly as follows: A key time-scale is the time-scale t c h ar for modifying the background distribution. Classically 
this time-scale is set by the inverse bounce-frequency, i.e. we have t 2 har ^ uj^ 2 oc <f> -1 . Mathematically this simply 
reflects the fact that the coupling to go and the higher harmonics g n is linear in the wave amplitude, as seen in Eqs 
When quantum effects enters the velocity derivative is replaced by a finite difference. Classically when 
the scale lengths in velocity space increases (cf Fig. 1) the magnitude of the nonlinear terms increases. However, 
the finite value of Sv q prevents the continuous increase of the nonlinear coupling, as is seen from Eqs <[T5])-([T7|), and 
effectively the transition from classical to quantum regime correspond to a change <f> <&/(! + Sv q ). Physically this 
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means that the characteristic time for acceleration of particles increases when the minimum velocity change comes in 
steps of 2 Sv q . Noting that 7 £^c har > 1 f° r the nonlinear modifications of go to occur faster than the linear damping 
explains the significance of the parameter R n \. 



V 
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Figure 5. The profile of go(v, t = 10) shown for initial amplitude <!>(£ = 0 ) = 4 for Sv q = 0.5 (lower panel) and for Sv q = 4 (upper 
panel). The discrete nature of the veloctiy dependence with peaks at discret values of Sv q is clear in the quantum regime. 
There is also a relation to the discrete energy states of Fig. 2, since the discrete velocity structure disappears when there is 
a large number of trapped states, as seen in the lower panel. However, there is not a one-to-one correspondence between the 
energy states and the velocity structure. Specifically the discrete structure of the velocity dependence can be present in the 
strong quantum regime when the particle trapping is absent due to quantum effects. 

The results reported here are similar in some respects to the findings of Ref. [ 21 ] and Ref. [ 22 ]. In Ref. m 
the full Wigner equation was solved numerically. It was then found that the nonlinear regime of Landau damping 
was suppressed when hk/mv t = 8 . This is a rather extreme quantum regime, and our condition for suppressing the 
nonlinear regime (essentially oo 2 mp being negative) which roughly gives hk 2 /2m > uj 2 b /^(l is much easier to fulfill. 
The nmain reason for our relaxed quantum condition is that we have focused on a resonance in the tail of the 
distribution, where 7 ^ u. In this case 5v q become comparable to the velocity scale length close to the resonance 
long before Sv q becomes comparable to the thermal (or Fermi) velocity. Thus we can investigate a regime where the 
wave-particle interaction is quantum mechanical, although the real part of the wave frequency is determined by the 
classical dispersion relation. By contrast, for a resonance in the bulk of the distribution, all quantum effects appear 
simultaneously. Next we compare our results with those of Ref. [22], that has studied the initial evolution of the wave 
particle interaction. We note that that the comparison of the quantum time-scales t q = 2m/hk 2 with the classical 
bounce time t B = 2tt/uj b made in [22 is very similar to the discussion made above. In particular Ref. [22] notes 
that the condition t q t B implies the absence of trapped particles. However, since Ref. [ 22 ] only studies the initial 
evolution, the conclusion that the nonlinear regime is suppressed for t q <t B suggested in that work is not accurate. 
Using the notation of [22], rather the condition for suppressing the nonlinear regime is t q < t 2 B ^L, as explained above. 

Much of the above features can be summarized in Fig. 7 that shows the different regimes of Landau damping 
plotted as a function of the nonlinearity parameter ujb/il and the quantum parameter hk 2 /2m^ B . The interesting 
part of the diagram is the nonlinear quantum regime u) B /jl > 1 and hk 2 /2m^ B >1. It should be noted that the 
condition for quantum suppression of the nonlinear amplitude oscillations given by Ref. [22] is the same condition 
that separates our region V from region VI (i.e. the condition that controls the existence of trapped particles). 
Interestingly, however, we find that the system can undergo nonlinear bounce-like oscillations even without trapped 
particles, and hence the region of quantum suppression (region IV) is determined by a distinct condition involving 
the quantum modified bounce frequency. It should be noted that the collisional influence can modify the different 
regimes to some extent. This will be discussed in some detail in the final section. 

A feature of Fig. 4 not explained by the above discussion is the small but continuous decrease of wave energy seen 
for the highest amplitude in Fig. 4. Besides the amplitude oscillation there is a continuous decrease in 4>(t) for initial 
condition 4>(t = 0) = 40 but not for <f>(t = 0) = 20. To understand this we must investigate how the energy of the 
particles in the resonance region is distributed among the harmonics. As described above the convergence is the sum 

over harmonics is fast for R tr = \f&/5v q <C 1. For the curves displayed in Fig 4 only the one with <F(£ = 0) = 40 has 
significant harmonic generation. For all the other curves the harmonics g n for n > 2 is effectively not excited, and 
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the two terms (1/2) f res gj)dv+ f res \g 2 \ dv accounts for more than 99% of the particle energy. However, due to the 
rapid phase mixing of g<± for large Sv q , this harmonics does not return the particle energy to wave energy, but instead 
the second harmonic acts as a leakage of energy, not taking part in bounce-like oscillations. This can be contrasted 
against the oscillatory evolution of the energy of the third harmonic for Sv q = 0.5 which could be seen in Fig. 6. 
To a small extent the continuous loss of wave energy to higher harmonics can be observed also for <F(£ = 0) = 20, 
but to a much smaller degree. This is illustrated in Fig 8, where the relative energy content in the second harmonic 
/res |#2 | dv/Wtot is plotted for <F(£ = 0) = 20 and for <!>(£ = 0) = 40. As can be seen, by a comparison with Fig. 4 
the continuous increase of energy in the second harmonic accounts rather well for the overall downward trend in wave 
energy. A similar mechanism in principle exists also for lower values of Sv q . However here it is much less effective. 
The reason is twofold. Firstly, for lower Sv q the peaks of #2 and higher harmonics occurs for a smaller velocity, and 
phase mixing is less effective. Thus the energy of the higher harmonics can be returned back to lower harmonics more 
easily. Secondly, for the same degree of nonlinearity (same value of R n \) the trapping parameter is larger for smaller 
Sv q . Hence a larger number of harmonics are excited for lower 5v q . This means that the highest harmonic excited 
takes a smaller proportion of the energy, in which case the leakage of energy is smaller. In practice, the details of the 
energy loss is likely to be determined by collisional effects, as very small scales in velocity space develops during the 
evolution [25]. 



Figure 6. Comparison of the relative energy in the third harmonic f \g 3 (t)\ 2 dv/Wtot for initial amplitude 4>(t = 0) = 4. The 
two curves are computed for 8v q = 0.5 and for 5v q = 4 respectively. 


V. SUMMARY AND DISCUSSION 

Usually particle dispersive effects as accounted for by the Wigner-Moyal equation are assumed to be significant 
when the thermal de Broglie length AdB = h/mv t , or its counterpart for a degenerate plasma, h/mv f (where up is the 
Fermi velocity), becomes comparable to the Debye length Ap = Vt/oo p (or up/coy). The physics behind this condition 
is that the importance of collective effects require scale lengths not much shorter than the Debye length, and for the 
sinus-operator in Eq. ^ not to reduce to its classical limit we need hL~^ ar /mv c har ~ 1- If the characteristic spatial 
scale length L c p ar is assumed not to be smaller than Ap, and the velocity scale length u c h ar is of the order of v t we 
get the condition hoj p /mv 2 ~ 1 (or hbj p /mv\ ^ 1 for a degenerate plasma) for quantum effects to be significant [l - J. 
This condition applies broadly to many situations, but does not hold for wave-particle interaction in general, as we 
have u c har vt in case the resonance lies in the tail of the distribution. In the present paper, focusing on the case 
of Langmuir waves of a single wavelength, we have found that quantum effects become important in the nonlinear 
regime once hk 2 /m ~ ujb in which case the trapped particle energies are quantized. Already when hk 2 /m > yp 
nonlinear oscillations at the bounce frequency is slowed down in accordance with ujb ^b/{ 1 + Sv q ) 1 ' 2 . While these 
oscillations are qualitatively similar to the well-known bounce oscillations |S9], it should be noted that in the quantum 
regime such oscillations can occur even in the absence of trapped particles. Increasing the quantum parameter 5v q 
even further, eventually the nonlinear oscillations are suppressed when ujb/{ 1 + dvq) 1 ^ 2 < yp in which case usual 
linear damping takes place. 
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Figure 7. Summary of the different regimes of Landau damping for a resonance in the tail. Region I and II: Linear Landau 
damping. The quantum effects in region II does not change the value of the linear damping for a resonance in the tail Region III: 
Classical nonlinear regim with bounce oscillation of frequency ~ ujb- Region IV: Quantum suppression of nonlinear bounce-like 
oscillations. Region V: Nonlinear regime with bounce-like oscillations in the absence of trapped particles Region VI: Region of 
trapped particles with quantum modified bounce oscillations. Note that the same expression ~ ujb/( 1 + hk 2 /tti^l) 1 ^ 2 for the 
quantum modified amplitude oscillations apply in both region V and VI 



Figure 8. Comparison of the relative energy in the second harmonic / \g2{t )\ 2 dv/Wtot for Sv q = 20. The two curves are 
computed for initial amplitude 4>(t = 0) = 20 and &(t = 0) =40 respectively. 


Since the phenomena of study evolves on a time scale slower than the plasma frequency, the above picture can be 
modified when the collisional influence is accounted for. Let us illustrate this by considering a concrete example. First 
we note that the long term evolution will always be affected by collisions. The characteristic time scale for collitions 
to be important r c is given by min[z/ e [ 1 , z/” 1 ]. To some extent we can assure that the phenomena of study occur on a 
faster scale by picking a relatively large amplitude and avoding a resonance too far out in the tail of the distribution. 
However, for a strongly coupled plasma this will not suffice, as z/ ee approaches the plasma frequency rather than being 
a small parameter. Thus for our model to apply we must first of all consider a plasma that is weakly coupled. As 
an example we pick a plasma with a number density no — 10 28 m -3 and a temperature T ~ 10 7 K that correspond 
to a coupling parameter T ~ 0.01. Such parameter values may result from laser-plasma experiments, see e.g. Ref. 
m- Furthermore the temperature is well above the Fermi temperature, and we have a Debye length of the order 
Xd ~ 10 - 9 m. Adjusting the wavenumber we can chose the amount of resonant particles that determine the linear 
Landau damping rate and thereby pick he quantum parameter. In order to get Sv q ^4 we can aim at 7 l — 10 12 s -1 
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which is obtained for a wave number k ~ 2 x 10 8 m _1 . As 5v q > 1 we are well in the quantum regime. We note, 
however, that although the plasma is weakly coupled, the time-scale for collisional damping is still shorter than the 
the linear damping time y^ 1 11 . This is typically the case in the quantum regime, and hence the long-time evolution 
is generally much affected by collisions. However, most of our findings concern the nonlinear quantum behavior. 
As we will demonstrate below these findings are to a large extent still relevant even when collisions are taken into 
account. Firstly, the quantum modified bounce frequency ujb —» cjb/(1 + fivq) 1 / 2 can be experimentally verified after 
a few bounce oscillations. For example, if we pick 4> > 10V the bounce oscillations take place on a much faster 
scale than collisional damping, in which case the quantum modification of the bounce frequency becomes detectable. 
Moreover, generally the prediction of a nonlinear regime in the absence of trapped particles (regime V in Fig. 7) can 
be verified after a few (quantum modified) bounce oscillations. Hence in an experimental situation we do not need 
to follow the evolution until collisional effects become important. The situation is somewhat different when it comes 
to the condition for quantum suppression of the nonlinear oscillations. It is possible to verify that previously given 
conditions are incorrect in agreement with our theory, as the presence of bounce oscillations at a faster scale than the 
linear damping rate will confirm this. However, the quantitative confirmation of our condition that separates regions 
IV and V in Fig 7 is likely not possible, as that would require us to study the system when the quantum modified 
bounce frequency is of the same order as the linear Landau damping rate. In this case we need to follow the evolution 
long enough such that collisions will modify the picture. Thus we conclude that most features of Fig 7 remains even 
when we account for collisional effects, but that there are restrictions arising from the collisional influence that makes 
the boundary between region IV and V uncertain. 

While the study here has focused on Landau damping due to electrons, it should be noted that a similar process 
may take place for photons ESI EH- Photon Landau damping can be relevant to Langmuir waves with relativistic 
phase velocities, in the presence of electromagnetic radiation, when resonant electrons are nearly absent. The wave 
nature of the photons provides the classical counterpart to the quantum electron states, and we can similarly identify 
a photon bounce frequency uJBph ^ wb- Photon trapping effects were experimentally observed by [32] . 
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